Petrological evidence for shock‐induced high‐P metamorphism in a gabbro          下载免费PDF全文

J.‐J. Yang  H.‐R. Zhang  A.‐P. Chen  M.‐X. Huang 《Journal of Metamorphic Geology》2017,35(2):121-140

Microlites (minute spherulitic, dendritic, skeletal, acicular and poikilitic crystals) diagnostic of crystallization in quenched melt or glass in fault rocks have been used to infer fossil earthquakes. High‐P microlites and crystallites are described here in a variably eclogitized gabbro, the wallrock to the coesite‐bearing eclogite breccia at Yangkou in the Chinese Su‐Lu high‐P metamorphic belt. The studied hand specimens are free of discernible shear deformation, although microfractures are not uncommon under the microscope. In the least eclogitized gabbro, the metagabbro, stellate growths of high‐P minerals on the relict igneous minerals are common. Dendritic garnet crystals (<1?5 μm) grew around rutile and/or phengite replacing ilmenite and biotite, respectively. Skeletal garnet also rims broken flakes of igneous biotite and mechanically twinned augite. Radial intergrowths of omphacite and quartz developed around relict igneous orthopyroxene and are rimmed by skeletal or poikilitic garnet where a Ti‐bearing mineral relict is present. Acicular epidote, kyanite and phengite crystallites are randomly distributed in a matrix of Na‐rich plagioclase, forming the pseudomorphs after igneous plagioclase. In the more eclogitized gabbro, the coronitic eclogite located closer to the eclogite breccia, all the igneous minerals broke down into high‐P assemblages. Thick coronas of poikilitic garnet grew between the pseudomorphs after igneous plagioclase and ferromagnesian minerals. The igneous plagioclase is replaced by omphacite crystallites, with minor amounts of phengite and kyanite. Thermodynamic modelling of the plagioclase pseudomorphs shows an increase in P–T in the wallrock from the metagabbro to the coronitic eclogite, and the P–T variation is unrelated to H₂O content. The fluid‐poor pressure overstepping scenario is unsupported both by phase diagram modelling and by whole‐rock chemical data, which show that the various types of eclogitized gabbro are all fairly dry. A large pressure difference of >2 GPa between the metagabbro and the coesite‐bearing eclogites ~20 m apart cannot be explained by the subduction hypothesis because this would require a depth difference of >60 km. The microlites and crystallites are evidence for dynamic crystallization due to rapid cooling because constitutional supercooling was unlikely for the plagioclase pseudomorphs. The lack of annealing of the broken biotite and augite overgrown by strain free skeletal garnet is consistent with a transient high‐P–T event at a low ambient temperature (<300 °C), probably in the crust. Therefore, the eclogitization of the wallrock to the eclogite breccia was also coseismic, as proposed earlier for the eclogite facies fault rocks. The outcrop‐scale P–T variation and the transient nature of the high‐P–T event are inconsistent with the other existing tectonic models for high‐P metamorphism. The fact that the less refractory but denser biotite is largely preserved while the more refractory but less dense plagioclase broke down completely into high‐P microlite assemblages in the metagabbro indicates a significant rise in pressure rather than temperature. Given that the metamorphic temperatures are far below the melting temperatures of most of the gabbroic minerals under fluid‐absent conditions, stress‐induced amorphization appears to be the more likely mechanism of the coseismic high‐P metamorphism.  相似文献   

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation          下载免费PDF全文

Chongbin Zhao  Thomas Poulet  Klaus Regenauer‐Lieb 《国际地质力学数值与分析法杂志》2015,39(6):571-593

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions          下载免费PDF全文

Xiong Zhang  Jean‐Louis Briaud 《国际地质力学数值与分析法杂志》2015,39(13):1369-1409

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Soil property variation mapping through data mining of soil category maps          下载免费PDF全文

Fei Du  A‐Xing Zhu  Lawrence Band  Jing Liu 《水文研究》2015,29(11):2491-2503

Spatial information on soil properties is an important input to hydrological models. In current hydrological modelling practices, soil property information is often derived from soil category maps by the linking method in which a representative soil property value is linked to each soil polygon. Limited by the area‐class nature of soil category maps, the derived soil property variation is discontinuous and less detailed than high resolution digital terrain or remote sensing data. This research proposed dmSoil, a data‐mining‐based approach to derive continuous and spatially detailed soil property information from soil category maps. First, the soil–environment relationships are extracted through data mining of a soil map. The similarity of the soil at each location to different soil types in the soil map is then estimated using the mined relationships. Prediction of soil property values at each location is made by combining the similarities of the soil at that location to different soil types and the representative soil property values of these soil types. The new approach was applied in the Raffelson Watershed and Pleasant Valley in the Driftless Area of Wisconsin, United States to map soil A horizon texture (in both areas) and depth to soil C horizon (in Pleasant Valley). The property maps from the dmSoil approach capture the spatial gradation and details of soil properties better than those from the linking method. The new approach also shows consistent accuracy improvement at validation points. In addition to the improved performances, the inputs for the dmSoil approach are easy to prepare, and the approach itself is simple to deploy. It provides an effective way to derive better soil property information from soil category maps for hydrological modelling. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Pre‐Alpine discordant granitic dikes in the metamorphic core of the Betic Cordillera: tectonic implications     

Antonio Sánchez‐Navas  Antonio García‐Casco  Agustín Martín‐Algarra 《地学学报》2014,26(6):477-486

Pegmatite dikes bearing andalusite crosscut foliation S₂ in Alpujarride gneisses and schists. Post‐S₂ andalusite is transposed by a foliation S₃, defined by fibrolite, which affects the dikes. The dikes represent highly differentiated granitic magmas with low REE and Zr contents and a positive Eu anomaly. U‐Pb SHRIMP dating of magmatic zircons provided Pan‐African ages (cores) and late Variscan ages (rims). However, U‐rich rims also provided metamorphic Alpine ages, supporting a polyorogenic tectonometamorphic history for pre‐Mesozoic Alpujarride rocks.  相似文献   

Early solar system aqueous activity: K isotope evidence from Allende     

Yun Jiang  Piers Koefoed  Olga Pravdivtseva  Heng Chen  Chun‐Hui Li  Fang Huang  Li‐Ping Qin  Jia Liu  Kun Wang 《Meteoritics & planetary science》2021,56(1):61-76

The alkali element K is moderately volatile and fluid mobile; thus, it can be influenced by both primary processes (evaporation and recondensation) in the solar nebula and secondary processes (thermal and aqueous alteration) in the parent body. Since these primary and secondary processes would induce different isotopic fractionations, K isotopes could become a potential tracer to distinguish them. Using recently developed methods with improved precision (0.05‰, 95% confidence interval), we systematically measured the K isotopic compositions and major/trace elemental compositions of chondritic components (18 chondrules, 3 CAIs, 2 matrices, and 5 bulks) in the carbonaceous chondrite fall Allende. Among all the components analyzed in this study, CAIs, which formed initially under high‐temperature conditions in the solar nebula and were dominated by nominally K‐free refractory minerals, have the highest K₂O content (average 0.53 wt%) and have K isotope compositions most enriched in heavy isotopes (δ⁴¹K: ?0.30 to ?0.25‰). Such an observation is consistent with previous petrologic studies that show CAIs in Allende have undergone alkali enrichment during metasomatism. In contrast, chondrules contain lower K₂O content (0.003–0.17 wt%) and generally lighter K isotope compositions (δ⁴¹K: ?0.87‰ to ?0.24‰). The matrix and bulks are nearly identical in K₂O content and K isotope compositions (0.02–0.05 wt%; δ⁴¹K: ?0.62 to ? 0.46‰), which are, as expected, right in the middle of CAIs and chondrules. This strongly indicates that most of the chondritic components of Allende suffered aqueous alteration and their K isotopic compositions are the ramification of Allende parent‐body processing instead of primary nebular signatures. Nevertheless, we propose the small K isotope fractionations observed (< 1‰) among Allende components are likely similar to the overall range of K isotopic fractionation that occurred in nebular environment. Furthermore, the K isotope compositions seen in the components of Allende in this study are consistent with MC‐ICP‐MS analyses of the components in ordinary chondrites, which also show an absence of large (10‰) isotope fractionations. This is not expected as evaporation experiments in nebular conditions suggest there should be large K isotopic fractionations. Nevertheless, possible nebular processes such as chondrules back exchanging with ambient gas when they formed could explain this lack of large K isotopic variation.  相似文献   

Does stress transmission in forelands depend on structural style? Distinctive stress magnitudes during Sevier thin‐skinned and Laramide thick‐skinned layer‐parallel shortening in the Bighorn Basin (USA) revealed by stylolite and calcite twinning paleopiezometry     

Nicolas Beaudoin  Olivier Lacombe  Marie‐Elonore David  Daniel Koehn 《地学学报》2020,32(3):225-233

The Sheep Mountain‐Little Sheep Mountain Anticlines, Bighorn Basin (USA) formed as basement‐cored Laramide structures in the formerly undeformed foreland of the thin‐skinned Sevier orogen. We take advantage of the well‐constrained microstructural network there to reconstruct differential stress magnitudes that prevailed during both Sevier and Laramide layer‐parallel shortening (LPS), before the onset of large‐scale folding. Differential stress magnitudes determined from tectonic stylolites are compared and combined to previous stress estimates from calcite twinning paleopiezometry in the same formations. During stress loading related to LPS, differential stress magnitudes transmitted from the distant Sevier thin‐skinned orogen into the sedimentary cover of the Bighorn basin (11–43 MPa) are higher than the differential stress magnitudes accompanying the early stage of LPS related to the thick‐skinned Laramide deformation (2–19 MPa). This study illustrates that the tectonic style of an orogen affects the transmission of early orogenic stress into the stable continental interior.  相似文献   

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock          下载免费PDF全文

J.‐W. Park  C. Park  J.‐W. Song  E.‐S. Park  J.‐J. Song 《国际地质力学数值与分析法杂志》2017,41(6):880-898

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

The influence of coupled physical swelling and chemical reactions on deformable geomaterials     

Yue Ma  Xiao‐Hui Chen  Lee J. Hosking  Hai‐Sui Yu  Hywel R. Thomas  Simon Norris 《国际地质力学数值与分析法杂志》2021,45(1):64-82

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.  相似文献   

Dating the onset of Lower Tagus River terrace formation using electron spin resonance          下载免费PDF全文

PIERLUIGI ROSINA  PIERRE VOINCHET  JEAN‐JACQUES BAHAIN  JORGE CRISTOVÃO  CHRISTOPHE FALGUÈRES 《第四纪科学杂志》2014,29(2):153-162

  相似文献